NACA RM 57B25 


6 / 



RESEARCH MEMORANDUM 


- v 


PROPAGATION OF SOUND INTO A WIND-CREATED SHADOW ZONE 

By David C. Pridmor e-Brown 
Massachusetts Institute of Technology 


library copy 

APR 29 1957 


LANGLEY AEBONHimCU UBORATOR ' 

LIBRARY NACA 

LANGLEY FIELD, VIRGINIA 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

April 22 , 1957 




1A 



NACA EM 57B25 


3 1176 01434 40 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


RESEARCH MEMORANDUM 


PROPAGATION OF SOUND INTO A WIND-CREATED SHADOW ZONE 
By David C. Pridmore -Brown 


SUMMARY 


The general wave equation is derived governing the propagation of 
sound in a stratified mo ving medium, the velocity of which varies only 
along one coordinate. Under the assumption that the flow velocity is 
small and slowly varying, a simplified equation is adopted which is sat- 
isfactory for the present application. A solution is found to this equa- 
tion corresponding to the pressure field around a spherical source located 
above a plane ground in a horizontal wind whose velocity increases with 
altitude. It is shown that within the acoustic shadow that forms on the 
upwind side of such a source the pressure field is similar to that which 
is obtained in the corresponding problem of a sound source in a tempera- 
ture gradient. An expression is derived for the rate of attenuation 
within a shadow which is brought about by the presence of both a wind and 
a temperature gradient. 


INTRODUCTION 


An important phenomenon associated with the propagation of sound in 
the atmosphere is that of the refraction of the acoustic rays which can 
be brought about not only by the presence of a temperature gradient but 
also by the presence of a wind with a gradient in speed. A uniform, wind 
will have little effect on the sound propagation since its speed will, in 
general, be much smaller than that of sound. However, the presence of a 
wind gradient leads to an effective variation in the speed of sound at 
different points, and this causes a refraction of the rays, just as in 
the case of a temperature gradient. The important difference between the 
two phenomena lies in the fact that in a temperature gradient the local 
speed of sound is determined by the temperature at that point, whereas in 
the wind case the effective sound speed depends on the direction as well 
as on the magnitude of the wind velocity. This means that, assuming the 
gradient in each cane to lie always in the same direction, call it the 
vertical, the sound field around a spherical source in a temperature 
gradient will be synmetrical about that vertical which passes through the 
source; in other words, it will depend on only two coordinates. In the 
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wind case, however, this symmetry clearly no longer exists, and the field 
will necessarily depend on three coordinates. 

This phenomenon leads to an interesting type of shadow formation in 
the wind case. If a spherical source is situated above ground in a hori- 
zontal wind whose velocity increases uniformly with height, then a shadow 
region into which no acoustic rays penetrate will form on the upwind side 
of the source but not on the downwind side. In fact, the distance to the 
shadow will be a minimum directly into the wind and will increase to 
infinity at right angles to it. This situation is illustrated in fig- 
ure 1. In the presence of a negative temperature gradient alone, the 
shadow distance is, of course, the same in all directions. 

A consideration of these facts suggested the tentative assumption 
which was made in a previous report (ref. l) on the sound field in a wind- 
created shadow. There it was assumed that an expression for the dif- 
fracted field within the wind-created shadow could be obtained from the 
corresponding expression for the temperature-created shadow by replacing 
the sonic velocity gradient by the wind velocity gradient times the cosine 
of the angle between the sound and wind directions . The results of the 
present study furnish the justification for this assumption. 

This investigation was conducted at the Massachusetts Institute of 
Technology Acoustics Laboratory under the sponsorship and with the finan- 
cial assistance of the National Advisory Committee for Aeronautics. 


SYMBOLS 


An 

a 

Ci,C 2 

c 

D 

h 


(1)/ -iit\ 

constant in ^A R e J = 0 

defined in equation (ll) 
contours 

sonic speed 

solution to equations ( 5 ) or (l8) which represents down- 
going wave 

Hanltel function of second kind of order n 

source height 

imaginary 
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J 0 ( x ) 

k 

M 

P 

Re 

r 


zero-order Bessel function 
propagation constant, cd/c 
Mach, number of wind, V/c 
sound pressure 
real 

radial coordinate parallel to ground plane, distance from 
source 

distance from source at height h to shadow boundary at 
height z 


s — dM/dz 


t 


U 


u 


V 

w 


x,y,z 1 

x l, x 2, x 3j 

z 

a 

P 

7 = 9-0 

6 = ±1 


time 

travel time from source at height h to shadow boundary 
at height z 

solution to equations (5) or‘(l8) which represents upgoing 
wave 

particle velocity; function of z, fc , and 7 defined in 
equation (22) 

wind velocity 

defined in equation (13} 

rectangular coordinates 

acoustic impedance of ground 

wave -number component in x-direction 

wave -number component in y-direction 




impedance ratio, z/p Q c 
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9 polar angle of wave number, tan“^(p/a) 

k wave -number component in horizontal plane, \Ja?- + 3 2 

p = 1 or 2 

p density fluctuations 

p Q density of medium 

0 polar angle 

cd angular frequency 


ANALYSIS 


Derivation of Wave Equation 

In an inviscid medium moving with the velocity V the linearized 
Navier-Stokes equations for the conservation of mass and momentum become, 
in the absence of sources. 


|p + V |£- + 

U"b uXj 



0 


(1) 


^ u i , .. ^( V5 il) , v Su i 1 dp 

dt 5 Bx-j cix-^ p Q ^x.j_ 


provided the wind velocity is assumed to lie in the x-^-direction and to 

be a function only of x^. Here p Q represents the static density of 

the medium which Is assumed constant; the density fluctuations p, the 
sound pressure p, and the particle velocity u^ are considered to be 

small quantities of the first order. If the further assumption is made 
that p = c 2 p, that Is, that the pressure and density variations are 
adiabatically related, then these equations yield 


^2 

_1_ ^_P 

c 2 St 2 


- M 2 ') 


*2 

d_p 

cbc 2 


a 2 p t a 2 p 

dy 2 dz 2 


2M 5 2 p 
c dx dt 


2 Po c 


dz dx 


( 3 ) 
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where M = V/c and the coordinates Xj_ have been relabeled x, y, 
and z. 


Xt is next assumed that the sound pressure p and the vertical 
component of* the particle velocity uz vary harmonically in time as 
-ikc-h 

e and can be Fourier analysed in the x- and y-directions : 


p = e -ikct JJ giCccx+py^^p^^ dj3 




= e -ikct JJ e i( ' aX+Py ^G(a^ J z)dct dp 

The Fourier transforms are related by equation ( 2 ) which gives 


(*0 


AT? \ 

= ip 0 c(k - oM)G 
dz 


Substitution of these relations into original differential equation (j) 
then shows that F must satisfy 


£s.+ _ 2 °MX 3E + Rk-aMja 

^ 7 2 k - qM dz L 


a 2 - B 2 ]f = 0 


( 5 ) 


together with appropriate boundary and source conditions . Here 
M' = dM/dz. 

It is convenient to rewrite expressions (k) in polar coordinates by 
putting 


x = r cos (f 
y = r cos 0 
a = k cos 0 
P = k sin '0 
da dp = k. dk dS 
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so that p becomes, with 7=0 - 0, 

p = JJ~ e iKr COS ^ F(k,7,z)k 6k d.7 (6) 


The boundary condition at the ground surface, z = 0, can be speci- 
fied by requiring that the ratio of the pressure p to the normal com- 
ponent of the particle velocity u in this plane be equal to a complex 
impedance Z independent of the angle of incidence . A point source of 
sound is considered to be located a height h above this plane. 

The source condition can be derived in the usual way by writing p 
as the sum of two functions, namely, p = p Q + pp, where p Q represents 

the solution in an unbounded medium and pp is an everywhere -regular 
solution to wave equation ( 5 ) chosen so that p Q + p-j_ satisfy the 
boundary conditions. Since one of the requirements on p Q is that it 

represent outgoing radiation at great heights on the upwind side of the 
source, it will be convenient to label the two independent solutions of 
equation ( 5 ) as D(z) and U(z), where D(z) represents a downcoming 
wave and U(z) an upgoing one for large values of z. One can then 
write p Q in the form given by equation (6), where the corresponding 

F Q is given below the source by 


F 0 = AD(k,7,z) 


(7a) 


and above the source by 


F 0 * BU(k,7,z) 


(7b) 


jz-h I 

which corresponds to e"^ ' in a stationary medium. The two 

constants A and B are determined by the conditions at the source, 
namely, by requiring that the pressure be continuous across the plane of 
the source, z = h, and that the particle velocity (or pressure gradient) 
suffer a discontinuity across this plane, which is determined by the 
source strength. The first of these conditions yields the relation 


AD(h) = BU(h) 


( 8 ) 



NACA EM 57B25 


7 


The second source condition may he obtained by specifying the source 
output. To do this requires that 



dS = 4at 


(9) 


■where u^ is the camponfent of u perpendicular to the surface ele- 
ment dS; the integration is carried out over a small surface enclosing 
the source. Now, from momentum equation ( 2 ) there follows 


-ip Q c(k - aM)uj 


PqCUjM'Sj! = 

J 


on making use of equations (4) . The three components of the particle 
velocity are then 


where 


U 1 = 




Uq = a 


u* = a 


Sxq 

$P 

Sxq 

dp 

5x^r 


5x7 


a” 1 - ip e(k - aM) 




b" 1 = - oM) 2 


Po C . 


( 10 ) 


( 11 ) 


Thus 


J ^ dS = Ja(' ; 7p) n dS + h dxg dx^ 
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The last Integral Is over the two surfaces X]_ = 0“ and x^ = 0 + , 
and It vanishes since p and Its derivatives are continuous along x^. 

However, because of the discontinuity in the vertical component of Vp, 
the first integral on the right does not vanish but becomes 





the integration being taken over the two surfaces x^ = h~ and xj = h + . 
Thus, introducing an expression of the form given by equation ( 6 ) for p Q 
and reverting to polar coordinates r, z gives 


ircr cos /^'(h) _ AD'(h}j k dx d? r dr 


which must be equated to Here equations ( 7 ) have been used to 

express F Q , and the primes denote differentiation with respect to z. 

From the well-known relations 




iter cos 
e 


7 = J 0 ( Kr ) 


J 0 (Kr)K dK = 6 (r) = 0 



dr = 1 


(r / 0 ) 


There follows immediately 


a { bU' (h) - AD'(h)] = i. 


( 12 ) 


where a^ = a(h) . 
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The two equations (8) and (12) determine A and B and, hence, 
also the unbounded solution p Q , which becomes above the source (z > h) 


Po = d* J '* dr e lltr COS 7 £ D(h)U(z) (15) 


Here 


W = lp Q c jk - ctM(h)] (U'D - D'U) 


(3A) 


and is independent of z. A s imi lar expression holds below the source 
but with z and h interchanged . 

The general solution p 0 + pp can now be written 


P 



iKr cos r iji(h) + F^U(z)k <3k dr 


(15) 


where Fp l 0 determined by requiring equation ( 15 ) to satisfy the normal 

V z 

impedance boundary condition — — - — at z = 0 . It turns out to 

grad p iop G 


be 


*i 


KD(O) - HD'CO) 
kU(O) - i£U'(0) 


U(h) 


( 16 ) 


where 
to z. 



and the primes again denote differentiation with respect 


Approximate Solutions Valid Near Ground 

Up to now no restrictions have been placed on the Mach number M 
beyond requiring that it vary in only one direction. In the general case 
approximate solutions to equation ( 5 ) can be readily obtained by a modi- 
fied W.K.B. method (cf., e.g., the work of I anger , ref. 2). For the 
present application it will be sufficient to assume that both and 

M’ are so small that the terms involving these quantities in equation (5) 
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can be neglected. It will also be assumed that M increases linearly 
with altitude from the value zero at the ground, z = 0, so that 


M = sz 


( 17 ) 


•where s is a constant. 

With these restrictions equation (5) becomes 


aSe + (kS . a 2 . g2 

dz 2 x 


2kasz 


) r - 


0 


or 


d 2 ? 

dz^ 


+ 



K 2 


- 2kKsz cos 



(18) 


The two restrictions on the range of validity of this equation can 
be expressed as 



s/k « 1 

( 19 ) 

and 

sz « 1 

(20) 

The first is a high-frequency requirement, while the second is a 
tion to low heights; in practical cases in atmospheric acoustics 
restriction seems to be very stringent. 

res trie 
neither 


The general solution to equation (l8) is of the form 



F(k, 7 ,z) = u u 3 / 2 ) 

(21) 

where 

-2/V 

u = |k 2 — k 2 - 2kicsz cos (7 + 0)J j2ekKs cos (7 + 0)J 

(22) 
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and represents a linear combination of one -third-order Bessel 

functions. The quantity e - ±1 -will be specified later. The func- 
tion U(z) can be written provisionally in terms of a Hankel function 
of the n-th kind as 


u(z) 



( 23 ) 


where = 1 or 2 depending on the value of e and will be specified 
later. The function D(z) is then defined similarly in terms of the 
other Hankel function. 


Transformation of Integral in Complex Plane 

It Is of interest to notice that if the function F(k,7, z) = [b(h) + 
F-jJ U(z ) in equation (15 ) did not depend on 7 then the integration over 

7 could be carried out directly and would yield J 0 (Kr), which was the 
starting point in the problem of a sound source In a temperature gradient 
(ref. J). In that case it proved advantageous to divide up the Integral 

Into the sum of two integrals by writing 2 J 0 = H^ 1 ) + . Guided by 

these considerations the 7 contour in equation ( 15 ) is deformed Into 
the sum of two Hankel function contours as shown in figure 2 . This Is 
allowed since the Integrand involves 7 only in the form cos 7 and 
cos (7 + 0 ), which are unchanged If 7 is Increased by 2 *, so that 
clearly the contributions from the infinite branches of the contours 
cancel each other. As in the case of the integral representation of the 
Hankel functions the integral converges provided that 


-T| < arg nr < jt - q 


and it. 

The 

solution 

^ 00 



O 

■tT 

iKr 

e 

cos 7 k , 
W 

00 



0 

& 

iKr 

e 

cos 7 jt 1 


W 


(210 
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The convergence “Of* these integrals in the complex k plane is next 
examined with a view to deforming the k integration path into closed 
contours, as was done in the temperature problem. In other words, for 
the k integration In the first Integral in equation (24), a contour 
might be taken enclosing the first quadrant in the clockwise direction, 
and, similarly, a contour enclosing the fourth quadrant in the counter- 
clockwise direction for the second integral as indicated in figure 3 - 
To justify this it must be shown that the contributions from the infinite 
arcs vanish, that there are no branch points in the first and fourth 
quadrants of the k plans, and that the Integrals along the imaginary 
axes cancel each other. 


The 

functions 
that is, 


absence of branch points follows directly from the fact that the 
a are integral transcendental functions of u. 




.3 ~ ) 

they are single valued and analytic for all finite values of 


u. 


(They are the 
finite values 


solutions to 
of u.) Also, 


(d 2 F/du 2 ) + uF - 0, which Is regular for all 
v u(k) is analytic in the first and fourth 


quadrants of the complex K plane, and, hence, e 1 * 1 " 003 7 F(k, 7 ,z) is 
free of branch points in this region. 


The behavior of the integrands on the infinite arcs Is next studied. 

First of all it is noticed that the exponential term e^ Kr cos 7 tends to 
zero as nr tends to infinity In the first quadrant for values of y 
lying on the contour and similarly this exponential tends to zero as 

Kr tends to infinity in the fourth quadrant for values of y lying on 
the contour Cq. It must also be shown that F(k, 7 ,z) -» 0 as Kr -»« 

when ReK > 0 (Re indicates real values). In the limit of large values 
of K, u becomes 


u 


K 2Pl 


1 + 2 ~ sz cos 




2 €kKS cos 


(y + 0 )] 

(r + 0)] 2 ^ 5 


whence 


|u5/2 


(*i)e 


3 ks cos (7 + 0 ) 


+ KZ 


If the following convention is adopted 

1/2 


jk 2 - k 2 - 2 krcsz cos (7 -f- 0)J = ijjt 2 - k 2 + 2 kKsz cos (7 + 0/j 


1/2 
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for k 2 > k 2 - 2kKsz cos (7 + 0 ), then. It is necessary to choose the minus 
sign. This requires that p = 1 if e = -1 and p = 2 if e = 1. It 
is then easily found from the asymptotic forms of equation ( 23 ) that 


D(h)U(z) 



(u3/2- ul 3/a) 


.1 if 2e Ss. 

* K 2 


1/3 


cos(7 + 0} 


-ic(z-h) 


where u^ = u(h) . Similarly, 


kD(O) - i£D'(0) 
kU(O) - I^U*(0) 


,-iei t/6 ( 


exp 




3ks cos (7 + 


U(z)U(h) 


I i 26 cos(7 + 0) 

* [_ K 2 


1/3 


e -2K 2 /3ks cos(7+0) e -K(z+h)^53ti€/l2 


whence 


f(k,7,z) 


i 1/5 r 

26 ^ cos (7 + 0 ) le -14 ^ 


z-h) + e -tc(zH*hjj 


and so has the required behavior. 

It can now be proved that the two Integrals taken along the imagi- 
nary axes cancel each other, that Is, that 


n 

dO J Cq 


iter cos 


^ F(k,7)k dK d 7 + /" F e lKr cos ^ F(k,7)k dK d7 = 0 
d 0 d C Q 


( 25 ) 
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for 


— T| < arg ter < it - q 

where q lies between 0 and it. To do this kt is replaced by 
rcre ijt in the second integral in equation (25) to obtain 



iKr cob 7 F( .^ r)K a* a r 


with the requirement -q < arg Kre ijt < it - q or -q - it < arg ter < -q. 
The Integration path C 2 ' (shown by the heavy dashed line in fig. 2 ) 

is now the same as Cq taken in the opposite sense and displaced by it. 
Accordingly 7 is replaced by 7 + it, so that this integral takes the 
form 



< 3 * d7 


with 


-q < arg kt < it - q 

For equation (25) to hold it is then necessary for 


F(ic, 7 ) = f[-k ,(7 + it J 

Reference to equation ( 22 ) shows that this condition is fulfilled in the 
present case and, hence, that equation (25) is valid. 

The general expression for the pressure field p has now been 
reduced to the sum of two contour integrals taken around the first and 
fourth quadrants of the k plane. These Integrals can be expressed as 
a sum of their residues taken at the poles of F(k,7,z) in this plane. 
Finally, the integrations over 7 can be evaluated approximately by a 
saddle-point integration through the points of stationary phase of the 
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Integrands; namely, 7 = 0 for the contour and 7 = k for the 

C2 contour . 

First of all, on carrying out the approximate integration over 7, 
equation (2 1) becomes 


J e 1 ^ ^^3 E F(k,0,z) dK + 
(1st quadrant) 

— e _i £ f(k,jt,z) dK 

\Jk w 

( 4 th quadrant) 




(26) 


At this point the analysis is restricted to the case of a perfectly 
hard ground by letting £ -» 00 in equation (16) . The more general case 
can be treated by the methods used in the temperature problem (ref. k ) . 
For a hard ground the poles of F are given by the roots of U' (0) = 0 , 
which becomes 

H -2/ 3 (f u3/2 ) - ° 


whence 

| u 0 3/2 = Jy, 1 * 


where the components are constants, = 0 . 686 , Ag = 3-90, 

Aj = 7.O5, • • • , A n = jn - (3/^)] « (where n is large). 

Substituting for Uq gives 

(k 2^ K n g)5 /a _ 2 £W 

2ekKs cos (7 + 0) 2 
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or 

2/3 

K n 2 = k 2 - I^SKjjSkAn cos (7 + 0) e £i ^j 

As a result of the restriction s/k « 1 (eq. ( 19 ) ) , this equation can 
he writen approximately 



It is now required that ImK n > 0 (where Ixn indicates imaginary 
values) for < 0 < (damped outgoing propagation on upwind side of 

source) . This leads to the result that e = -1 and p. = 1. 

Since the values of K n all lie in the first quadrant, the contour 

Integral enclosing the fourth quadrant must vanish; this leaves the 
integral enclosing the first quadrant in the k plane, which can be 
evaluated as a sum of residues in the form 


P 



[*nr-(*A)],j( Ka/Z ) 

* n D '( K n'°) 

W / "N \ 


(f-V (*n,0) 
\OK/ 


(28) 


By using the Wronsklan relationship D'U - DU* = u^W, which gives 
D' = a^W/U at k = rc n , equation ( 28 ) can be rewritten in terms of U 

only as follows 


P 


V l 8jtK n [ a-LU(h)U(z) 

L V r 6 U(0)(d/dK)U'(0) 

n=l L 


(29) 


t '- 1/ V 1) (! u3/2 ) 


where 
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= ip 0 c(k - KSh cos 0) 


High-Frequency Behavior of Solution 

As in the corresponding temperature problem, the high-frequency 
"behavior of the solution (eq. (29)) can be studied in the shadow zone by 
replacing the functions U by their asymptotic forms, which are valid 
for large values of u. It is then found, on reintroducing the time 

factor e" ic0 k, that 


= ^B n r" 1 / 2 exp li l^r + + - ootjl 

n=l ^ J 


(32) 


where 


B - l 2kk 




H l/3 


( 1 ) 




-2 


On substituting for from equation (31) into equation (30) 

u = (-2| cos ^ 1 / 5 kz + A^e -1 ’ 1 ) 2 ^ 5 

is obtained which gives approximately, for large values of kz. 


2u 5/2 ^ + 25 n (z/a)^-/^j 
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where 


a = - 2 s cob f 


8 n * §('3§ ^ ) 2/5 e _Si ”/ 3 


then 


The terms in the exponent of the asymptotic solution (?2) become 


V + f[^G 3/2 + f &(h)]5/ 2 - tub = K n r + k 1^- o^/ 2 ( B 5/a + h 5/2) + 

SSuCz/o) 1 / 2 + (h/a) 1 / 2 - rntj 
= k(r - r D ) (1 - 6 n) - c(t - t Q ) 


Here, these quantities have been identified with r Q , the horizontal 
distance from source to shadow boundary, and with t Q , the corresponding 
travel along the limiting ray, which have been derived from ray acoustics 
in the appendix. 


In terms of these quantities the high-frequency form of the solu- 
tion is 


p (Constant) n 


rt -1 


(s/k) 1 / 6 -1/2 -^(t'to-h^l-Sn)! 

w 1 e 


(33) 


Within the shadow zone the solution is adequately expressed by the 
first mode alone because of the rapid attenuation of the higher modes. 
To this approximation the sound pressure is damped at the rate of 

8 . 68 lmk 62 6 .l(-s cos ^)^/^k^^db/unit distance (34) 

within the shadow ^ < f < r > r G j . This is in addition to the 

damping due to cylindrical divergence, which is expressed by the 
factor r - ^^ . 
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The combined effects of a wind and a temperature gradient can be 
accounted for by allowing k to vary with z. In the high-frequency 
approximation this results in the factor -s cos 0 in equation (34) 


being replaced by -s cos 0 - — — , where — is positive for a tem- 

c dz dz 

perature inversion. A shadow will form in this caBe within that sector 


for which the factor 


^-s cos 0 - i is real. By replacing 


by i , the result is obtained that a shadow will form within the 
c dz 

sector defined by 



provided that 0 is now measured from the 180° line ( -x-axis ) and 


*£> o. 

dz 


If 


dc 

dz 


dV 

dz 


the velocity gradient. 


then the temperature gradient predominates over 
In this case there will be a shadow in all 


directions if 


dc 

dz 


< 0 and no shadow at al 3 if 


> 0 . 

dz 


Within the normal or illuminated zone the high-frequency (ray 
acoustics) approximation to the solution can be readily obtained by 
integrating equation (24) around the points of stationary phase of the 
integrand in the k plane (ref. 5). 


Massachusetts Institute of Technology, 
Cambridge, Mass., April 6, 1956. 
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APPENDIX A 

RAY ACOUSTICS IN A WIND GRADIENT 


The 

requires 


acoustic -ray paths are determined by Fermat's principle, which 
that they be such that the integral 


/ 


ds 

c + V sin 0 


cos 0 


(Al) 


is stationary. Here ds is an increment of path length along the ray, 
V is the wind speed, 0 is the angle of inclination of the ray from 
the vertical, and 0 is the constant polar angle between the ray and 
the wind direction. 

From the geometry of the problem it is seen that 


sin 6 = r 

/l + (r ' ) 2 


where r' = With this substitution, and writing ds = \ll + (r') 2 dz, 

CLZ ? 

the Euler equation is obtained for which equation (Al) is stationary in 
the form 


or 



1 + (r 1 ) 2 ' 

1 + (r 1 ) 2 + Mr' cos 0, 


= 0 


— . 1 t (r llg = K (A2) 

Jl + (r 1 ) 2 + Mr' cos 0 


where M = V/c and K is a ray parameter. The limiting ray which 
defines the shadow boundary can be specified by requiring that r' -» *> 
at the ground, where the wind velocity is zero. This requirement gives 
K = 1 for the limiting ray. Carrying out the differentiating in equa- 
tion (A2) and solving for r' gives 



HACA BM 57B25 


21 


r« = 1 (l + S M cos 0 + . . (A3) 

y -2M cos 0 ' ^ / 

if hi ghe r powers of M are ignored. If M is set equal to sz it is 
found to the first order in M that 


r 



dz 

\j- 2sz cos 0 




■2s cos 0 


The total horizontal distance between the source at height h and a 
point at height z on the shadow boundary is then 



(A4) 


The corresponding travel time t Q of the limiting ray between 
these two points is given by the integral in equation (Al), which becomes, 
on substituting for ds and sin 0 in terms of r’, 

kf , x i (r ' )2 «■ 

yl + (r 1 ) 2 + Mr' cos 0 


Equation (A3) is now used to write r 1 = 1/V-2M cos 0, and M is set 
equal to sz. Then, to the same approximation as before, 

cto = r Q + j(-2s cos 0) 1 / 2 ^ z 3/2 + h.2/2^ (A5) 


is obtained, where r Q is given In equation (A^) . 
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Figure 1.- Illustration of mechanism of shadow-zone formation. 
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